Wave function mapping of graphene quantum dots with soft confinement 
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Using low-temperature scanning tunneling spectroscopy, we map the local density of states 
(LDOS) of graphene quantum dots supported on Ir(lll). Due to a band gap in the projected 
Ir band structure around the graphene K point, the electronic properties of the QDs are domi- 
nantly graphene-like. Indeed, we compare the results favorably with tight binding calculations on 
the honeycomb lattice based on parameters derived from density functional theory We find that 
the interaction with the substrate near the edge of the island gradually opens a gap in the Dirac 
cone, which implies soft- wall confinement. Interestingly, this confinement results in highly symmet- 
ric wave functions. Further influences of the substrate are given by the known moire potential and 
a 10 % penetration of an Ir surface resonance into the graphene layer. 

PACS numbers: 73.20.At, 72.10.Fk, 73.21. Fg, 73.22.Pr 



Graphene has moved in short time from first prepa- 
ration as a small flake [1] towards possible applications 
such as high frequency transistors [2] , supercapacitors [3] 
or touch screens [4]. Another exciting perspective is to 
use graphene quantum dots (QDs) as spin qubits [5]. The 
basic prerequisite is a very long spin coherence time [6], 
which might exist in graphene [7] due to the absence of 
hyperfine coupling in isotopically pure material and the 
small spin-orbit coupling [8]. First graphene QDs have 
been produced and probed by transport measurements 
[9, 10]. However, since graphene provides no natural gap, 
it is difficult to control the electron number [11]. More- 
over, the 2D sublattice symmetry makes the QD prop- 
erties very susceptible to the atomic edge configuration 
[5] unlike conventional QDs. As a result, chaotic Dirac 
billiards have been predicted [12] and were even claimed 
to be realized [9, 13], i.e. the wave functions are assumed 
to be rather disordered. To achieve improved control of 
graphene QDs, the QD edges must be well defined and 
a deeper understanding of the QD properties is manda- 
tory. 

Direct insight into QD properties is provided by scan- 
ning tunneling spectroscopy (STS) which maps out the 
squared wave functions of QDs [14] and, at the same 
time, determines the shape of the QD atom by atom. Us- 
ing STS, we probe graphene QDs with well defined zig- 
zag edges supported on an Ir(lll) surface [15]. These 
QDs maintain graphene properties as the filled part of 
the graphene Dirac cone lies in the Ir projected band 
gap [16]. By comparing the measured wave functions 
with model calculations, we determine the relationship 
between geometry and electronic properties and extract 
general trends. Most notably, the soft edge potential pro- 
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vided by the interaction of the QD edges with the sub- 
strate enhances the geometrical symmetry of the wave 
functions, thus rendering the QD more regular. The sus- 
ceptibility of the wave functions to the edge configuration 
is intimately related to the additional sublattice symme- 
try (pseudospin) which makes graphene so special [17]. 
Also the moire pattern induced by the graphene-Ir lat- 
tice mismatch [16] and the hybridization of graphene with 
an Ir surface resonance are shown to have an influence 
on the measured wave functions. 

STM measurements are performed in ultrahigh vac- 
uum at T = 5K [18]. Monolayer graphene islands are 
prepared by exposing clean Ir(lll) for 4 min to a pres- 
sure of 10 ~^ Pa of C2H4 at 300 K and subsequent an- 
nealing to 1320 K (30 s) [19]. The resulting graphene 
QDs have diameters of 2 — 40 nm as shown in Fig. la. 




FIG. 1. (color online) (a) (100 x 100) nm^ STM image of 
Ir(lll) covered by monolayer graphene islands; U = —0.3 V, 
/ = 0.3 nA; (b) atomically resolved (12 x 12) nm^ image 
of graphene island; (c) magnified view of zigzag edge with 
graphene lattice overlaid; U — 0.7 V, / = 20 nA. 
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Atomically resolved QD images (Fig. lb— c) reveal the 
complete enclosure of the QDs by zigzag edges. 

The local density of states (LDOS) of 15 islands 
is mapped by STS. We use a lock-in technique with 
modulation frequency v = 1.4 kHz and amplitude 
^ mod = IQniV resulting in an energy resolution ^ 
V'(3.3 • ksTy^ (1.8 • eUmod)^ = 18meV [20]. For 
dl/dU curves, we stabilize the tip at sample voltage 
t/gtab and current /gtab- Figure 2a shows a dl/dU curve 
laterally averaged over the hexagonal QD shown to the 
right. It displays three maxima below the Dirac point 
£^D, which is slightly above the Fermi level Ep [16]. 
Thus, the peaks belong to confined hole states. Fig. 2b- 
d show dl/dU maps at the peak energies. For the first 
peak {U = —0.26 V), one maximum of the LDOS in the 
center of the island appears, a ring shaped structure is 
observed at = —0.42 V, and, a maximum-minimum- 
maximum sequence from the center towards the rim with 
an additional star-shaped angular dependence is visible 
at = —0.63 V. We checked that no other LDOS shapes 
are present at —lAeV <U < V. From the sequence of 
observed LDOS shapes we conclude that they represent 
confined states of the QD. 

To model the QD states, we employ third-nearest neigh- 
bor tight binding (TB) calculations [21-23] using the 
atomic configuration of the QD found by STM, 

^ = Y1 1^^'^) ^(^'i) 1^^'^) ((/>j,s| + ^.c.. (1) 

i,s {ij),s 

The 7(i,j) are hopping amplitudes between sites i and j 
being = (3.14,0.042,0.35) eV for the (first, second, 
third) nearest-neighbors [21]. The Vi represent local on- 
site potentials. 

We first employed a spatially constant Vi within the is- 
lands, i.e. hard- wall-confinement. Regular, but also very 
irregular wave functions result, as shown in Fig. 2h and 
Fig. 3e-g. The irregular wave functions often display a 
large intensity at the rim of the QDs and illustrate the 
sensitivity of graphene QDs to details of the edge configu- 
ration [5, 12]. Such irregular shapes, however, were never 
found in the present STS experiments featuring about 50 
different states [24]. 

This failure is related to the two experimental facts that 
(i) a graphene fiake bends downward from D = 3.4 A in 
the center of a QD to D = 1.6 A at its rim [15, 26] and 
that (ii) the entire graphene fiake features a moire type 
corrugation leading to minigaps [16, 27, 28]. 

To incorporate effect (i) we determined the band struc- 
ture of graphene by ab initio density functional theory 
(DFT) calculations [24, 25] for different graphene-Ir sur- 
face distances D. Upper and lower limits for D were set 
by the known distance between extended graphene lay- 
ers and Ir(lll), D = 3.4 A [26] and the smallest distance 
found at the edge of a graphene island, D = 1.6 A [15]. A 
proper description of Ir(lll) surface states requires thick 
slabs which makes it unfeasible to use the large 10 x 10 
supercell necessary to account for the graphene-Ir lattice 
mismatch. Therefore, a slightly compressed Ir lattice is 
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FIG. 2. (color online) (a) black line: dI/dU{U) curve spa- 
tially averaged over the graphene QD shown to the right; 
t/stab = 0.5 V, /stab = 0.5 nA, [/mod = 10 mV; grey line: 
DOS(E) of the same island as obtained by TB calculation 
(see text); vertical bars mark the calculated eigenstate ener- 
gies with degeneracies indicated as numbers; (b)-(d) dl/dU 
images recorded at energies E = U ■ e as marked; / = 0.2 nA; 
Umod = 10 mV. (e)-(g) LDOS maps calculated with soft edge 
potential at energies indicated; (h),(i): LDOS of an individual 
state calculated without (h) and with (i) soft edge potential. 



used making graphene and Ir(lll) commensurate. This 
allows us to work with a slab of 24 Ir layers with graphene 
on both sides and a vacuum space of 20 A between slabs. 
The insets in Fig. 3a exhibit the resulting band struc- 
tures for two different fixed D. The size of the gap AEd 
is plotted in Fig. 3a. We incorporate the effect of the D 
dependent band-gap on Vi within the TB through [12]: 



Vi. 



AE[D{n)]/2-a, 



(2) 



where the Pauli matrix acts on the sublattice degree 
of freedom. A homogeneous Vi^rim would open a gap 
of size AE at £^d- The functional form of AE[D] = 
(0.7* (3.6- L)[A])2 + 0.23) eV is taken from the fit to the 
DFT calculations (Fig. 3a). We model the global height 
variation of a graphene QD by linear increase of D{r) 
from the rim towards 10 A inside the island as suggested 
by the DFT calculations of [15]. We checked that reason- 
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FIG. 3. (color online) (a) Energy gap versus graphene-Ir 
distance D as deduced from DFT calculations; insets: band 
structure around for two different D as marked by arrows 
with indicated; grey area: projected bulk bands of Ir; 
Thick black lines: graphene states; (b)-(g) Calculated LDOS 
(= 1^1^) for individual confined states with energies marked: 
(b)-(d) with soft edge potential; (e)-(g) without soft edge po- 
tential; (h) experimental kn • r — En/{hvB) ■ r for the two 
peaks closest to Ed at different average island radius r; cir- 
cles: n = 0, squares: n = 1; dotted lines: zeros of the first 
two Bessel functions (see text). 



able modifications do not change the results significantly 
[24]. 

To incorporate effect (ii), we added a moire potential Vi^^n 
to Vi. Based on the experimentally observed minigap of 
200 meV [16, 27, 28], we use a harmonic variation of m 
in each of the three dense packed directions of graphene 
with a total amplitude of 400 meV [24] . Finally, the peak 
width r of the eigenstates is adapted to the experiment 
leading to V{E) = 0.33 • \E\. 

The resulting LDOS curve (grey line. Fig. 2a) as well 
as the calculated LDOS maps (Fig. 2e-g) exhibit excel- 
lent agreement with the experimental data. Importantly, 
the calculations yield only states that reflect the hexag- 
onal symmetry of the QD shape in agreement with ex- 
periment, but none of the irregular states found without 
smooth confinement [24] . This can be rationalized by the 
suppressed interaction of the confined states with the zig- 
zag edges, which would break sublattice symmetry [29]. 
The increased geometrical symmetry is illustrated in Fig. 
3b-g comparing wave functions of the same quantum dot 
with soft (hard) confinement leading to symmetric (ir- 
regular) states. Thus, softly opening a band gap at the 
QD edge leads to strongly improved control on the states 
residing in its interior. 

To illustrate this crucial finding, we show that the state 
energies in our QDs can be correctly estimated by a sim- 
plified circular flake geometry. We obtain En = hvj^^kn 
with Dirac velocity = 10^ m/s and kn deduced from 



the Bessel functions: 

Jn(kn-r) = 0, n = 0,l,.... (3) 

Up to an island area oi A = 150 nm^ (average radius: 
r = y^/tt), the estimate fits the experimental peak 
energies to within ~ 20 % for the two lowest energy 
states (Fig. 3h). Larger islands do not follow this trend 
because of their strong deviation from a circular shape 
(e.g. Fig. 4a). Obviously, neither the sensitive sublat- 
tice symmetry of graphene [5], nor the influence of the 
iridium substrate enter Eq. (3) showing the simplicity of 
softly confined graphene QDs. Note, in addition, that 
the agreement in Fig. 3h only uses the peak energies as 
an experimental reference and does not refer to the mea- 
sured LDOS shapes. Thus, peak energies are compatible 
with = 1 ±0.1 • 10^ m/s. 

In larger islands, we observe the influence of Vi^^ on 
wave function patterns directly, at energies E < — 0.6eV. 
Figure 4a shows an STM topography of a large QD ex- 
hibiting a regular moire pattern [19]. The dl/dU map 
in Fig. 4b and the calculated LDOS in Fig. 4c repro- 
duce the moire topography albeit with inverted ampli- 
tude. The same result is found for all larger islands [24]. 
We checked that normalizing the dl/dU images to ac- 
count for a spatially varying tip-surface distance [30] did 
not change the LDOS patterns. 

One feature, already visible by comparing Fig. 4b and c, 
is not accounted for by a spatially varying Vf. a bright 
rim of the island in the dl/dU image. This rim is found 
for all islands, but cannot be reproduced by the TB cal- 
culations [24]. Closer to £^d, this feature develops into 
a standing wave pattern that finds its counterpart out- 
side the island with slightly larger wave length A (Fig. 
4d— f). The dispersion relations E{Ak = tt/A) [31] inside 
and outside the islands are evaluated as displayed in Fig. 
4e and in [24] for 11 islands. They are shown together 
with results from standing waves at step edges of Ir(lll) 
in Fig. 4g. The E{Ak) curves are linear according to 
E -hv^Ak + Eb with ^ 4.9 • 10^ m/s, Eb -0.3 
eV outside the island and — 4.5 • 10^ m/s, E^ = —0.2 
eV inside the island. These values agree with those of the 
Ir surface resonance So around F found by photoemission 
(dashed line) including the energy offset between the two 
E{Ak) curves [28]. The values disagree with vb for the 
graphene Dirac cone on Ir(lll) by a factor of two and 
with Er> for the Ir S2 surface state by 0.5 eV [16]. Thus, 
the standing wave patterns within the QD are attributed 
to an intrusion of So into graphene. The amplitude of the 
standing wave in the islands Aq is found to be close to the 
amplitude outside the island Ai^ for several islands and 
energies [24]. This is surprising considering the fact that 
the tip is 0.23 nm further away from the Ir surface, when 
positioned above graphene, which would suggest a reduc- 
tion in dl/dU intensity by a factor of 100 [20]. However, 
DFT calculations reveal that So, exhibiting sp-symmetry, 
penetrates into graphene. The ratio between the LDOS 
in the graphene layer Iq and the LDOS in the Ir surface 
layer /jr is i^DFT = Ic/hr S - 12 % (inset of Fig. 4h). 
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FIG. 4. (color online) (a) STM image and (b) dl/dU map 
of a large graphene QD; 30 x 30 nm^ U = -0.65 V, / = 
0.5 nA, [/mod = 10 mV; (c) calculated LDOS of the same 
QD at = -0.65 eV; (d)-(f) dl/dU maps of a graphene QD 
recorded at the energies marked; 27 x 30 nm^, / = 0.5 nA, 
Umod = 10 mV; deduced wave lengths Aout (Ain) outside (in- 
side) the QD are marked in (e); (g) resulting dispersion rela- 
tions E{Ak = vr/Ain/out) inside (stars) and outside (triangles) 
of the QD as well as from standing waves scattered at Ir(lll) 
step edges (circles); full lines are linear fits with resulting 
v-D indicated; energy offset is marked; dashed line is deduced 
from photoemission on clean Ir(lll) [28]; (h) relative inten- 
sity R of So and S2 in graphene as deduced from STS data 
(squares) and from DFT calculations (So: circles, S2: trian- 
gles); inset: calculated LDOS of So at ^ = —0.4 eV along the 
direction perpendicular to the surface; /ir and Ic as used for 
determination of R are marked. 



For comparison, S2 shows only Rbft — 0.02 %. Figure 
4h favorably compares i^DFT of So with the data from 
STS i^sTS where the apparent Aq/^if is rescaled accord- 
ing to i^sTS = Ac/Air ■ e"^ [20] with a = 1.1 - 1.2/ A 
deduced from I{z) curves and S = 1.1 A being the differ- 
ence between real height (3.4 A [26]) and apparent STM 
height (2.3 A) of the graphene above the Ir(lll). Thus, 



we can quantitatively reproduce the strength of So intru- 
sion into graphene. A simple explanation for the strong 
So intrusion is not obvious, but we note that, according to 
DFT, also the d^-like surface state Si, located at Ep and 
exhibiting no dispersion [16], penetrates into graphene 
with c:^ 10 — 40 % and the 7r-electrons of graphene pen- 
etrate back into Ir with R 1 — 4 %. 
Finally, we would like to comment on the fact that the 
So state partly dominates the LDOS patterns, while the 
peak energies are reproduced nicely by the Dirac cone of 
graphene. We assume that the life time of the graphene 
states is large enough to lead to confinement resonances 
appearing as peaks, while the life time of the So is sig- 
nificantly shorter leading only to exponentially decaying 
standing waves at the step edges of the graphene islands. 
Indeed, we do not observe peaks within the spectroscopy 
of the islands, where the standing wave of So is dominat- 
ing the LDOS pattern. Moreover, the standing wave gets 
always significantly weaker in intensity away from the 
step edge. Of course. So probably influences the LDOS 
patterns of the small islands as well, which might ex- 
plain the remaining deviations between theory and exper- 
iment in Fig. 2 (e.g. (d) and (g)). This subtle interplay 
between graphene electrons and So electrons within the 
graphene island might also explain the too low ^'d = 6-10^ 
m/s resulting from the analysis of the LDOS pattern of 
graphene quantum dot states on Ir(lll) in ref. [32]. 

In conclusion, we mapped the LDOS of graphene QDs 
supported on Ir(lll). For small islands, properties of an 
isolated graphene QD with soft edge potential reproduce 
the measured wave functions. Most importantly, the soft 
edge induced by the substrate is required for the experi- 
mentally observed high symmetry of the wave functions. 
Larger islands show an additional standing wave pattern 
caused by an intruding Ir surface resonance and signa- 
tures of the moire potential. 

We acknowledge helpful discussions with N. Atodire- 
sei, C. Stampfer, G. Burkard, S. Runte, and a referee, 
as well as financial support by DFG (LI 1050/2-1, MO 
858/8-2, BU 2197/2-1), Fonds National de la Recherche 
(Luxembourg), and FWF (SFB-F41 VICOM). Numeri- 
cal calculations are performed on the Vienna Scientific 
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Note added in proof: During the referee process, two 
publications with similar experimental results have been 
published [32], which were submitted later than our 
manuscript. 
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I. IMAGE GALLERY OF INVESTIGATED QUANTUM DOTS 



Figure 1 shows an STM image gallery of all investigated graphene quantum dots (QDs) 
on the Ir(lll) substrate. Detailed comparisons of dl/dU images and dl/dU curves with 
the results from third nearest neighbor tight binding calculations have been performed for 
theses islands, but only representative data are shown. 




FIG. 1. Constant current STM images of investigated graphene QDs. Tunneling parameters are 
marked in the images. 
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II. DENSITY FUNCTIONAL THEORY (DFT) CALCULATIONS 



Ab initio DFT simulations were carried out using the plane- wave PWSCF code included 
in the QUANTUM-ESPRESSO package [1]. 

We firstly studied the relationship between the size of the band gap in graphene and the 
graphene-Ir surface distance ranging from D — 3.4 A to D = 1.6 A. A 1 x 1 cell with 
the graphene cell parameter was used, which corresponds to a slightly compressed Ir(lll) 
surface. We use a slab of 24 Ir layers with graphene on both sides and a vacuum space of 
20 A between slabs. We employed gradient-corrected exchange correlation functionals [2] 
and fuUy-relativistic ultrasoft pseudopotentials including spin-orbit interactions [3]. The 
wave functions were expanded in plane waves with a kinetic energy cutoff of 30 Ry and a 
charge-density cutoff of 300 Ry. 20 x 20 x 1 Monkhorst-Pack meshes [4] of /c-points were 
used for the integration over the Brillouin zone. 

The electronic band structure of the system was computed along the T, T" and S lines 
of the surface Brillouin zone, corresponding to the path f — ^ — M — fin the reciprocal 
space. The graphene states were identified by projecting the wavefunctions of the slab on 
the atomic wavefunctions centered on the C atoms with a threshold of 50 %. This allowed 
us to calculate the gap AE between the lower and upper graphene cone as a function of 
the graphene-Ir distance (see main text). The absolute value of the Dirac point energy with 
respect to the Fermi level at equilibrium distance {D — 3.4 A) slightly deviates (by 200 
meV) from photoemission experiments [5] . The band gap has also been estimated from the 
photoemission data to be 100 meV [5], i.e. 100 meV smaller that the one found by DFT at 
the equilibrium distance. We assume that these discrepancies are caused by the compression 
of the Ir(lll) within the DFT calculations, but does not affect the general trend of increasing 
the band gap with shortened graphene-Ir distance. 

Very recently, Varykhalov et. al [6] detected a large Rashba effect on a surface state 
of Ir(lll) (denoted as So in the main text) near the T point: the properties of this state 
were found to be hardly affected when the surface is covered with graphene. We accordingly 
carried out DFT simulations of both the clean Ir(lll) surface with a lattice constant of 
ao — 2.758 A, corresponding to the value of a/cc = 3.90 A for bulk fee Ir obtained from 
DFT calculations, and the graphene-Ir system with D = 3.4 A. For these calculations, 
we also used a 24-layer Ir slab and identified the Ir(lll) surface states by projecting them 
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onto the atomic orbitals of the surface and subsurface Ir atoms requiring more than 25 % 
of their weight located at these atoms. We found the surface state So, which is rather a 
surface resonance, in both cases. The band structure of the graphene-Ir system along the 
T — M direction is shown in Fig. 2 with the surface resonance marked. The penetration 
of the surface resonance So (as well as of the surface state S2) into the graphene layer was 
calculated by relating the weight on the C atoms Ic to the weight on the Ir surface atoms 
/ir using i?DFT = Iq/ hr (scc main text). 




FIG. 2. Band structure of the graphene-Ir system along the V — M direction. The dispersion of 
the So surface resonance is indicated by a red line. The state has been identified by requiring that 
the sum of the squares of the projections of the state on the orbitals of the atoms of the surface 
and subsurface layer be larger than 0.25. 

III. CONFINED POTENTIAL USED FOR THE CALCULATION 

The potential Vi used for the tight binding calculation is modeled according to the distance 
D from the graphene sheet to the Ir surface [see Eq. (2) in main text]. We use a Berry- 
Mondragon like potential Vi = AE[D]/2 • to model the gap opening at the Dirac point, 
as taken from our DFT calculations (see Fig. 3a of main text). Consequently, the potential 
conserves K-K^ symmetry [7] and models the (partial) hybridization of the carbon orbitals 
with the Ir states known from our DFT calculations. Since the graphene quantum dot 
approaches the Ir surface at the edges, we obtain a smooth edge confinement that is key to 
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the formation of states with high symmetry observed in the experimental LDOS patterns. 
Near the edges, we assume a hnear approach (as the most simple model) of the graphene 
sheet towards the Ir substrate from the equilibrium distance of 3.4 A to 1.6 A at the edge 
over a distance of 10 A. To assert that this choice does not influence our conclusions, we 
have performed calculations for several diflFerent functional forms for distance between the 
quantum dot and the Iridium surface keeping the outmost distance of 1.6 A fixed. Excluding 
unphysical, vertical kinks in the shape of the graphene flake, we flnd no noticeable changes 
in the wavefunction patterns (see Fig. 3). The variations in calculated resonance energies for 
different types of edge potentials are below 15 meV, which is smaller than the experimental 
energy resolution. 

Secondly, the lattice mismatch between Ir and graphene is taken into account by a spatially 
varying moire potential Vm{x^y)^ which consists of a suitable superposition of sinusoidal 
functions in accordance with ref. [8, 9] (see main text). In order to obtain the amplitude 
Knax of we perform tight binding calculations of the bandstructure of an inflnitely 
extended graphene sheet in the presence of for different Vmax, and extract the size of the 
minigap A induced by zone folding [5]. Comparing to ARPES measurements, which feature 
A ^ 200 meV [5], we obtain 14iax = 400 meV. 
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FIG. 3. The three lowest-lying eigenstates with their eigenenergies marked for different 
parametrizations of the distance development between the graphene quantum dot and the irid- 
ium substrate at the edge of the quantum dot: (a)-(c) linear interpolation oc x over 15 A, i.e. 
D{r±) = 3.4 - 0.12 • (15 - r±) • 0{15 - r±) with r±, Dink and 0{x) being the step function; (d)-(f) 
over 10 A, i.e. D{r±) = 3.4 - 0.057 • (10 - r^f/'^O • (10 - r^); (g)-(i) quadratic interpolation 
(oc x^) over 10 A, i.e. D{r±) = 3.4 — 0.018 • (10 — r±)'^ • 6^(10 — r±)] in the main manuscript, a linear 
interpolation over 10 A was used (see main manuscript, Fig. 2). 

IV. CALCULATED STATES WITHOUT SOFT EDGE POTENTIAL 

Figure 4i shows an STM image of a particular island with the corresponding dl /dU maps 
displayed in Fig. 4 e— h. We display all patterns of dl/dU maps found in the voltage range 
between U = —0.2 V and U = —0.6 V. For comparison, wave functions are calculated by the 
third nearest neighbor tight-binding calculation neglecting the soft edge, i.e. Vi = oo outside 
the island. Inside the island, the moire potential is maintained. Superpositions of squared 
wave functions at the energies corresponding to the dl /dU images are shown in Fig. 4 a— d. 
Only the two very similar states belonging to K and are superposed in Fig. 4a— c. In 
Figure 4d already 6 states, which have a very similar energy, had to be superposed. The 
shapes of the wave functions and the LDOS correspond to the measured dl/dU maps and 
the energies reasonably well, but the calculated shapes are less regular and more extended 
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FIG. 4. (a)-(c) calculated squared wave functions of the island shown in (i) at energies marked; 

(d) calculated LDOS of the same island (consisting of 6 wave functions) at the energy marked; 

(e) -(h) dl/dU maps of the island shown in (i) at different energies E — U - e d^s marked; / = 0.3 
nA, U^od = 10 mV; (i) STM image of the graphene QD; U = -250 mV, / = 0.3 nA; (j)-(k) 
calculated squared wave functions of the QD shown in (i) not found in experiment; all calculations 
are with abrupt edge. 



towards the edge. Moreover, the calculated wave functions displayed in Fig. 4 j— k do not 
resemble the experimental LDOS maps. These wave functions feature a strong weight at 
the edge of the QD. They are suppressed by the soft confinement caused by hybridization 
of the graphene Pz orbitals with the iridium substrate. 
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V. STATES DOMINATED BY THE MOIRE POTENTIAL 

For energies below -0.6 eV, the dl/dU maps are dominated by the moire potential (see, 
e.g., the three islands in Fig. 5). The experimental data (left images) correspond to the 
calculated LDOS maps (right images) being prone to the moire potential with amplitude 
400 meV. Neglecting this moire potential leads to more uniform wave patterns in the LDOS 
maps (not shown). Notice that a bright rim is visible in all experimental LDOS maps of 
Fig. 5 as well as in Fig. 4(f)-(h), which is not reproduced by the calculation and attributed 
to the Ir surface resonance around F (see main text). 



-0.74 eV 




-0.77 eV 


•. . . 33x33 nm2 

V 

• • 

-0.67 eV / 


-0.7 eV 



38x38 nm2 



-0.67 eV 



-0.7 eV 



FIG. 5. dl/dU maps (left images) of different graphene QDs imaged at energies below -0.6 eV 
(/ = 0.4 nA, ?7mod — 10 mV). The images are dominated by the moire potential and are largely 
reproduced by the calculation (right images). 
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VI. STANDING WAVES OUTSIDE THE GRAPHENE ISLAND 



Figure 6 shows dl/dU maps of three different quantum dots where the contrast is tuned in 
order to see the standing wave outside the QD. Obviously the wavelengths of these standing 
waves Aout is decreasing with decreasing energy. By line scans, as shown in Fig. 6m for 
the dl/dU map in Fig. 6j, we deduced the £'(7r/Aout) dispersion shown in Fig. 4g of the 
main text. Since the measured standing wave pattern exhibits half the wave length of the 
impinging Bloch wave, Ak — 7r/Aout is used such that the dispersion E{Ak) can directly be 
compared with data from angular resolved photoelectron spectroscopy (ARPES) [6]. 
The standing waves are also observed at step edges of the Iridium(lll) surface not covered 
by graphene. Figure 7a shows an STM image of the uncovered Ir(lll) surface with two step 
edges. The two white dots are larger adsorbates on the surface. The dl/dU images exhibit 
standing waves at the step edges whose wave length decreases with decreasing energy. In 
addition, remaining oxygen adsorbates on the surface are visible as black dots which induce 
an additional complicated scattering pattern on the terraces. The wave lengths Air of the 
standing waves at the step edges are determined by line scans averaging along the step edge 
and the resulting E^-k/Xi^) is also plotted in Fig. 4g of the main text. The symbols in Fig. 
4g of the main text exhibit a very similar steepness of E{ti/\x) {x = in, out, Ir) for all three 
measurements, but the absolute values are lowest for the pure Ir(lll) surface, slightly higher 
(about 50 meV) for the standing waves around the graphene QDs and the highest for the 
standing waves inside the graphene QDs being another 100 meV higher. The photoemission 
data [6] show the same energy shift of the dispersion of the state So by about 150 meV 
between uncovered Ir(lll) and Ir(lll) completely covered with graphene. 



9 




2 4 6 8 10 12 14 

lateral displacement (nm) 



FIG. 6. (a)-(l) dl/dU maps of three different graphene QDs recorded at the energies E — U - e 
marked; ?7mod = 10 mV; (a)-(d),(i)-(l) / = 0.2 nA, (e)-(h) / = 0.5 nA; all images exhibit a standing 
wave around the island; wave length decreases with decreasing energy; blue line in (j) marks the 
line for the line scan shown in (m); (m) line scan along the line shown in (j); the deduced wave 
lengths of the standing waves inside (Ain) and outside (Aout) the QD are marked. 
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FIG. 7. (a) STM image of the Ir(lll) surface with two step edges; U = -0.2 V; / = 0.3 nA; 
(b)-(d) dl/dU maps of the same area recorded at the energies E = U • e marked (/ = 0.3 nA, 
Umod — 10 mV); standing waves are visible at the upper and the lower side of the step edges. 

VII. ESTIMATE OF CONFINED ENERGIES 

As described in the main text, we estimate the energy of confined states by the zeros of 
the first two Bessel functions according to: 

Jn{kn'r)=0, n = 0,l. (1) 

with eigenenergies En = hv^km a Dirac velocity of graphene vd = l - 10^ m/s and Radius 
r of the island. The first two peak energies Eq and Ei of the experiment are determined with 
respect to the Dirac point E^ for different islands. Since the shapes of larger islands become 
more irregular (see Fig. 1), only the nine smallest islands are considered. The average island 
radius r is deduced from the island area Ahy r = a^/AJtt. The resulting /c^ • r is plotted as 
a function of r for the two peak energies closest to Ejj in Fig. 3 (h) of the main text. For a 
few of the islands, the energy of the 2^^ resonance in the dl/dU curve is not well defined due 
to broadening of lineshapes by finite state lifetime, and thus not considered in the present 
analysis. For the smallest island, only the first resonance energy lies within the Ir projected 
band gap. Obviously, reasonable agreement of the model with the experiment is found. The 
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discrepancy increases with radius which is attributed to the more non-circular shape of the 
larger islands. 
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